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We revisit Taylor dispersion in oscillatory flows at zero Reynolds number, giving 
an alternative method of calculating the Taylor dispersivity that is easier to use 
with computer algebra packages to obtain exact expressions. We consider the effect 
of out-of-phase oscillatory shear and Poiscuillc flow, and show that the resulting 
Taylor dispersivity is independent of the phase difference. We also determine exact 
expressions for several examples of oscillatory power-law fluid flows. 
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1. Introduction 

The standard Taylor-dispersion problem (Taylor (1953)) is to determine the ap- 
proximate long-time distribution of some passive tracer with density <fi in a channel 
or pipe Q = M x S, where S is some connected smooth one- or two-dimensional 
manifold. The governing equation in PDE form is 

^ + V-(V^--DV0)=O, (1.1) 

where t is time, V is the velocity field, and D is the diffusivity, usually with zero- 
flux boundary conditions on the walls. Here we focus on the case where S = [0, a] 
and the component of the velocity along the channel is a periodic function of time. 

The early work on oscillatory flows focused on flow in pipes (Aris (I960)), and 
much later work on oceanographic applications in unbounded flows (Young et. al. 
(1982)), but more recently there has been interest in biomechanics, notably Schmidt 
et. al. (2005), and the references therein, who studied consequences of oscillatory 
flows on nutrient transport in bone during physical activity and ultrasonic therapy. 
In these biomechanical applications, the flow is probably better modelled by a 
combination of shear- and pressure-driven flows, and natural questions are whether 
these contributions interact, and whether the phase difference between the two 
components is important. 

Also in this study we find an expression for the Taylor dispersivity that is 
more appropriate for use with computer algebra packages. This extends the work 
of Jansons and Rogers (1995), by providing a more accessible version of the solution 
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in a form that could be easily applied by people not familiar with the theory of 
stochastic processes, although it uses stochastic calculus in its derivation. 

Using computer algebra, we then find the general expression for the Taylor 
dispersivity at zero Reynolds number for sinusoidal shear and Poiscuillc flow with 
arbitrary phase difference, and show in particular that the Taylor dispersivity is 
independent of this phase difference. This extends the work of Claes and Van de 
Broeck (1991). We also find exact expressions for the Taylor dispersivity at zero 
Reynolds number for several pressure-driven power-law fluid flows. 

2. Dispersion in an oscillatory flow 

We suppose that the vertical diffusion Y is confined to the interval [0, a], which in 
ltd form is governed by 

dY t = adB t + v(Y t )dt + dL° t -dL a t , (2.1) 

where B is a standard Brownian motion (that is, it has diffusivity |), and L° and 
L a are local-time terms to keep the particle in the interval. (The local times change 
only when Y hits a boundary.) The above equation describes, in fluid mechanical 
terms, a particle with diffusivity D = cr 2 /2 and vertical drift v. For simplicity, we 
take a to be constant, as it is in almost all applications, though this condition could 
easily be relaxed, and we also, for simplicity, assume that v is bounded. 

As is usual in Taylor-dispersion problems, we neglect the effect of particle dif- 
fusion in the x direction, though it would not be difficult to include this effect if it 
were considered useful. Thus the horizontal displacement is given by 

X t = f u{s,Y s )ds, (2.2) 
Jo 

where u(t,y) is a periodic function of time, with period 2tt/uj say, and "=" means 
equal by definition. 

From (2.1), we see that the stationary vertical density q is given by 

q(y)dy ee km P[Y t G dy] = a cxp ( jT ^-dr^j dy, (2.3) 

where a is a normalization constant (so J Q a q(y)dy = 1). 

The part of the horizontal displacement due to the cross-channel mean of hori- 
zontal advection is trivial, and unimportant in the long-time limit as it is bounded, 
so we remove this contribution from the problem by replacing u(t, y) with u'(t, y) = 

u (t, y) - Jo q(v) u (^ y) d v- 

To find the variance of X t , apply Ito's formula: 

#(*. Y t ) = f y {odB t + dLl - dLl) + (Vg + dt, (2.4) 

for some function / to be identified. Now take |^(t, 0) = |^(t, a) =0 for all t to 
remove the local-time terms, and 

i 2 d2 f , d f , 9 f <(. \ to ,s 

^ W +V dy- + M= U( ' y) ' ( 5) 
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where / is taken to be an oscillatory function of t, with period 2ir/ui. This gives 
f(t, Y t ) - /(0, Y ) - J* yt(s, Y s )dB s + X' t , (2.6) 

where 

X' t = f u'(s,Y s )ds. (2.7) 
Jo 

The integral term can be rewritten, namely 

j'^ { s,y, )d B, = W^(f y( s,Y, } )\y (2.8) 

where W is a standard Brownian motion. Then the ergodic theorem tell us that in 
the limit t — > oo, 



First observe that 



lim E[X?]/t= lim E[Xl]/t, (2.10) 



and second that / is bounded. Thus, (2.6), together with the other observations, 
implies that the Taylor dispersivity is given by 



/>2rc / to p 

V= ]imE[X*]/t=- / / 



Furthermore, in the limit t — > oo, 



^^N(0,V) (2.12) 



in distribution. 



3. Without vertical advection 

Throughout this section, we assume that the vertical advection v = 0. Furthermore, 
suppose that the horizontal advection has the form 

u(t,y) =u 1 (t,y)+u 2 (t,y), (3.1) 

where U\ and u 2 are both periodic functions of time with period 27r/w, and, for 
all t and y, that u[(t,y) = u[(t,a — y), and u' 2 (t,y) = —u' 2 {t,a — y), where again 
the dashes denote that the y means have been removed. In this case, since (2.5) is 
linear, / can be split into /i and f 2 forced respectively by u[ and u' 2 , and inheriting 
the same symmetries in y. Thus (2.11) implies that 

V = V x +V 2l (3.2) 
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where T>\ and V 2 are the Taylor dispersivities respectively for u\ and u 2 , since 
q = 1/a and 

J ^* = - (3 - 3) 

as the integrand is an odd function of y — a/2. 

Note that (3.2) does not hold in general for arbitrary v, unless the periods of 
the two oscillations arc different. 



4. Special cases 

We now consider several example flows in which it is possible to find the Taylor 
dispcrsivity exactly. We begin by considering the way in which shear- and pressure- 
driven flows interact in the case of a Newtonian fluid, and then consider several 
cases of power-law fluid flows. In these examples, the vertical advection velocity 
v is zero. The former, in effect, provide solutions to the Hclc-Shaw cell model of 
the corresponding porous material system, namely a microscopic flow driven by 
shearing the whole porous material, and a pressure-driven flow for a fixed porous 
material. For biomechanical applications, see Schmidt et. al. (2005). 

In general we need to solve (2.5), and the algebra is very messy even in simple 
cases. If v, the vertical advection velocity, is non-zero, even if it is constant, the 
corresponding expressions for Taylor dispersivity V are too complicated to show 
here, even for simple horizontal flows. However, using symbolic algebra packages, 
it is possible to find simple exact expressions for V in some cases. 

For simplicity in the rest of this section, we choose units that make o = 1, and 
a = 1. Also in these special cases v = 0, so (2.5) reduces to 

W + l = »'('•")• <«' 

with §£(t,0) = §£(M) = for all t. Since v = and a = 1, we find q = 1, 
which slightly simplifies (2.11). With this non-dimensionalization, the time-scale 
for the cross-channel relaxation to occur is 1 , and the time-scale of the oscillations 
is In the following examples, u is proportional to cos(wi). Thus, when doing 
the computer algebra, we can consider u and / to be the real parts of functions 
proportional to exp(iwi). 

(a) u(t, y) = cos(ujt)y and v = 

For u(t,y) — cos(wt)y, using computer algebra, we find from (4.1) and (2.11) 

_ ^cos(i^) + ^cosh(^) — sin(^) — sinh(^) 
1 = 2i/ 5 (cos(i/)+cosh(i/)) ' ( ' ' 

where v = uj^ (see Figure 1). In the limit of small u), 
_ 1 31lu 2 5461w 4 

V h 1 . (4.3) 

60 45360 194594400 v ' 

In the limit of large oj, 
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Figure 1. Plots of V(uj)/V{0). The upper curve is for u(t,y) — cos(cot)y(l — y) and the 
lower curve is for u(t, y) = cos(tot)y. 

This is in agreement with the result of Claes and Van de Broeck (1991). This is 
almost, but not exactly, the result previously found by Young ct. al. (1982) and 
extended by Jansons and Rogers (1995) for the case of unbounded flow. It differs 
by a factor of 2 from the result of Young et. al. (1982), which is explained by Jansons 
and Rogers (1995), and lacks the oscillatory term of Jansons and Rogers (1995). 
The large uu result here is, as expected, the average of the result for the unbounded 
case. 



(b) u(t,y) = cos(ujt)y(l — y) and v = 

Similarly for u(t, y) = cos(uit)y(l — y), and solving for V using computer algebra, 
we find 



^ _ vcos(v) — ^cosh(^) — 3sin(^) + 3sinh(^) 
6^ 5 (cos(^) — cosh(^)) 

(see Figure 1). In the limit of small ui, 



(4.5) 



V 



1 



3780 1496880 583783200 



+ ■ 



(4.6) 



In the limit of large to, 



(4.7) 



This agrees with the result obtained by Claes and Van de Broeck (1991). 
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(c) A linear combination of the flows in the previous examples 

We now extend the work of Clacs and Van de Broeck (1991), and consider a 
flow driven by both a moving upper boundary and a pressure gradient, but not 
necessarily in phase, namely 

u(t, y) = Ui cos(ujt)y + U 2 cos(ut + ip)y(l - y), (4.8) 

where ip is a constant. From (3.2), (4.2) and (4.5), we find 

V = U?d 1 +U$d 2 . (4.9) 

Notice that this result does not depend on the phase shift, which could be important 
in some applications. 

(d) Power-law fluids 

The same approach can be used to quickly find exact expressions for some 
non-Newtonian flows. Here we consider power-law fluid flows, driven by pressure 
gradients. However, most of the results are too complex to give here, and are essen- 
tially useful only in computer programs. However, with the procedure given here, 
the reader can quickly recalculate them with the help of computer algebra. These 
flows are of the following form: 

u n {t, y) = cos(wt) (const. — \y — \\ n ), (4-10) 

where the constant is fixed by the boundary conditions which may be no-slip or 
partial slip. (Note that n is not the defining parameter of a power-law fluid, but it 
is more convenient here.) However the results for the Taylor dispersivity does not 
depend on the constant. 

Exact expressions for the Taylor dispersivity can be found for some values of n. 
However, the expressions for integer n seem to be by far the simplest, and we give 
a few examples of these below. We were unable to find the solution for general n. 

As these flows are symmetrical about y = | , we solve in a half range which means 
we avoid the singularity at y = i, which is problematic for computer algebra. 

In the following results, we do not take account of the relaxation times for the 
fluid, so we must assume that the time-scale of the oscillation is long compared 
with these. 

The case n = 1 is a rescaling of the Newtonian shear-flow result above: 
= v cos( |) + v cosh(|) - 2 sin(|) - 2 sinh(|) 

2^5 ( cos( | ) + cosh(f) ) • ) 

The case n = 2 is for a Newtonian fluid, already considered above (see (4.5)): 

V = d 2 . (4.12) 

The case n — 3: 

V = (9 iv (-80 + i/ 4 ) cos(i/) - v (-80 + v A ) cosh{v) 



+ 80 (-4 + v 2 ) cosh(^) sin(^) - 5 (-32 + 8 v 2 + v A ) sm(u) 

+ 80 (4 + v 2 ) cos(^) sinh(^) + 5 (-32 - 8 v 2 + v A ) sinh(i/))) 
/(160^ 9 (cos(^) - cosh(i/))). 



(4.13) 
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Figure 2. Plots of V(to) /V(0) for ti„ with n increasing from 1 at the bottom to 6 at the 

top. 



The case n = 4: 

V={v (-672 + z/) cos(» - v (-672 + z/) cosh(^) 

- 7 (-144 + 24 v 2 + v A ) sin(v) + 7 (-144 - 24 v 2 + z/ 4 ) sinh(z/)) (4.14) 
/(56 z/ 9 (cos(v) -cosh(i/))). 

The case n = 5: 

D = (5 (k (414720 - 13824 z/ + 5 z/) cos(^) 
+ j/ (-414720 + 13824 z/ 4 - 5 z/ 8 ) cosh(i/) 

- 17280 (-96 + 24 v 2 + v 4 ) cosh(^) sin(^) 

- 45 (18432 - 4608 v 2 - 768 z/ 4 + 48 z/ 6 + i^ 8 ) sin(z/) (4.15) 

+ 17280 (-96 - 24 z/ 2 + z/ 4 ) cos(^) sinh(^) 

+ 45 (18432 + 4608 v 2 - 768 z/ 4 - 48 z/ 6 + z/ 8 ) sinh(z/))) 
/(4608z/ 13 (cos(i/) - cosh(z/))). 
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The case n = 6: 

V=(9(v (11827200 - 54560 + 7 v 8 ) cos(^) 
+ i/ (-11827200 + 54560 v A - 7 v 6 ) cosh(^) 

- 77 (230400 - 38400 v 2 - 2560 v A + 80 u 6 + u 8 ) sm{v) (4.16) 
+ 77 (230400 + 38400 v 2 - 2560 v A - 80 u 6 + z/ 8 ) sinh(z/))) 
/(39424^ 13 (cos(z/) - cosh(i/))). 

Notice that in each case ui 2 V tends to a constant as w — > oo. 

There is little point in going beyond n = 6 as it is almost plug flow, with almost 
all the dispersion occurring near the boundaries. 

The above results suggest that in the case of integer n there is probably some not 
too complex generating function for the Taylor dispersivities, but MATHEMATICA 
was unable to work through the necessary algebra. However, since it would need 
computer algebra to unpack the generating function, one is probably better off 
computing the required Taylor dispersivity directly. 

A plot of all the results of this section are given in Figure 2. Notice that even 
though we have considered only integer n, the curves are sufficiently close that 
interpolation would be expected to well approximate intermediate values over the 
interesting part of the u> range. Also, given that the case n = 2 is known, the fact 
that the other plots are close to this, and occur in the expected order, is a good 
consistency check on the other results. 



(e) Comment on large and small lo limits 
From (4.1) we can see that, for small u), the approximate balance is 



dy 2 



2j6i = u '^y)> ( 4 - 17 ) 



and therefore the Taylor dispersivity is approximately the time-average of the non- 
oscillatory Taylor dispersivities for the instantaneous flows. In the sinusoidal case 
this is just i of the non-oscillatory Taylor dispersivity for the peak u. 
For large u>, the approximate balance is 

^ = «'(*. y). ( 4 - 18 ) 

giving an asymptotic Taylor dispersivity 

1 f 1 I dU \ 2 
V ~7T^ [-r(y)) d V> asw^oo, (4.19) 



2uj 2 J V dy 
where U is defined by u(t, y) — cos(uit)U (y) . 



5. Conclusions 

We have shown that the alternative method for determining the Taylor dispersiv- 
ity for an oscillatory flow lends itself to exact evaluation using computer algebra. 
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Figure 1 shows that for the two simple flows considered, namely an oscillatory 
simple shear flow and an oscillatory Poiscuille flow, both have a strong frequency- 
dependent Taylor dispersivity. These systems therefore could be used for measur- 
ing particle diffusivities more accurately by choosing a frequency that moves the 
corresponding Taylor dispersivities into the sensitive range. Alternatively, the same 
physical system could be used for separating particles of different diffusivities. Ex- 
perimentally the oscillatory simple shear flow is particularly easy to set up in a 
laboratory. 

One interesting direction for future research is applications to transport in 
bones, and other connective tissues, subjected to oscillatory stresses (for example 
due to exercise or during ultra sound therapy). It appears that physical processes of 
the type discussed here are essential for normal function, and a better understand- 
ing could result in more effective treatment. These applications however involve 
much more complex geometries, although would be expected to be qualitatively 
similar. The solutions here are in effect those for the Hclc-Shaw cell model of a 
porous material. Given we have shown that the Taylor dispersivities for shear- and 
pressure-driven flows do not interact in a channel flow, it is a natural question to 
ask if this carries over (at least approximately) to these biomechanical applications, 
where both effects are likely to be present. 

It is clear that many more non-Newtonian examples are possible with this ap- 
proach. However, it is less clear if there are any nice exact expressions for flows with 
a non-zero vertical (i.e. cross-channel) component. The intermediate expressions in 
these cases are much more complex, with more branch cuts, but there does appear 
to be hope. This is one area in which computer algebra is currently weak. 
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